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Abstract 

A priori error estimates are derived for /ip-versions of the finite element method 
for discontinuous Galerkin approximations of a model class of linear, scalar, first-order 
hyperbolic conservation laws. These estimates are derived in a mesh-dependent norm 
in which the coefficients depend upon both the local mesh size h K and a number p K 
which can be identified with the spectral order of the local approximations over each 
element. The results generalize those of Johnson and Pitkaranta to hp- methods. 


1 Introduction 

The discontinuous Galerkin method has received renewed interest as a higher-order scheme 
for approximating solutions to hyperbolic conservation laws. The method can be interpreted 
as a natural higher-order extension of finite volume methods while overcoming some of the 
difficulties associated with standard Galerkin or spectral methods. Originally studied for 
linear problems with constant coefficients and fixed-order approximations by Lesaint and 
Raviart [1] and for linear problems with variable coefficients by Johnson and his collaborators 


[2, 3], the method can be regarded as an elementwise application of the standard Galerkin 
method in which jumps on element boundaries are admitted naturally in the formulation. 
Since the usual condition of continuity of the solution along inter-element boundaries is not 
enforced, the element equations are decoupled, thereby eliminating the need for solving large 
systems of algebraic equations, even as the polynomial degree of the approximation increases. 
The solution in neighboring elements is mildly coupled through the flux across element 
boundaries. These element boundary fluxes are approximated using a numerical flux function 
which incorporates the hyperbolic character of the conservation law in much the same way 
as the finite volume method. The significant difference between the discontinuous Galerkin 
method and higher-order finite volume methods is that the coefficients in the polynomial 
approximation of the solution within an element are obtained by solving the conservation 
law and not by some post-processing of solution mean values. 

As with any higher-order method, special treatment is required to prevent oscillations in 
solutions which contain steep gradients or discontinuities. In the work of Cockburn, IIou, and 
Shu [4], the discontinuous Galerkin method was shown to be TVB (Total Variation Bounded) 
in the solution mean values, provided that the solution satisfied certain conditions. These 
conditions were used to construct a projection strategy for controlling oscillations in the 
mean values. Unfortunately, these conditions are not sufficient to eliminate oscillations in 
the pointwise values of the solution. Extensions of the projection ideas to the entire solution 
can be found in Bey and Oden [5] and Flaherty et. al. [6]. The numerical experiments of 
Bey and Oden [5] showed that with the discontinuous formulation, oscillations are confined 
to elements containing the discontinuity (in contrast to global oscillations resulting from the 
standard Galerkin or spectral methods) and that the order of accuracy of the method is p + 1 
in smooth regions when using uniform meshes with polynomial approximations of degree p. 

The discontinuous Galerkin method is ideally suited to adaptive hp strategies and to 
parallel computing. An a priori estimate of the error in the solution provides a basis for an 
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adaptive strategy since one then knows how the error behaves as a function of the mesh size 
and the degree of the polynomial approximation in an element. In this note, we derive an a 
priori error estimate for an hp version of the discontinuous Galerkin method. This estimate 
extends the previous work of Johnson and Pitkaranta [2] who analyzed an h version of the 
method with a fixed polynomial degree. 

2 Model Problem 

For simplicity, we consider a convex polygonal domain fl. The domain boundary dfi. with 
an outward unit normal vector n(;r) consists of two parts: an inflow boundary T_ to be 
defined below and an ’’outflow” boundary T+ = dfl \ T_. We consider the following linear 


scalar hyperbolic model problem, 

up + ait = f in fl C 7Z 2 (1) 

(3 ■ n u — (3 ■ n g on T_ (2) 

where / € T 2 (fl), g £ Z 2 (T_), j3 = (/?i,/? 2 ) T is a constant unit vector, up = (3 ■ Vu, 
a = a(x) is a bounded measurable function on f 1 such that 0 < ao < <z(x), and T_ = {x € 


df!|/3-n(x) < 0}. Note that while solutions to (1) may be discontinuous across characteristic 
lines x(s) defined by | ^ = (3, the solution is continuous in the direction parallel to (3. 

For / € Z/ 2 (fl), the space of admissible functions for solutions to (1) is given by V(Ct) = {u € 
Z 2 (fl ) \vp 6 L 2 (Cl)}. We note that the trace of functions in V(Cl) exist only in the direction 
/3; therefore, we further restrict fl so that f3 ■ n ^ 0 on dSl. If u is a solution to (1) with 
boundary values satisfying (2), then u also satisfies the following variational equality: 

f (up + au)vdx + / uv\f3 ■ n\ds = f fv + f gv\(3 ■ n\ds Vu € F(fl) (3) 

Jn J r_ J n J r_ 
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3 Notation and Preliminaries 


For a domain D in 7 Z 2 , let (v,iv)d = f D vwdx and ||v|||} = (v,v)d- Let || • || m> £) denote the 
norm in the usual Sobolev space 

The starting point for the discontinuous Galerkin method is (3) defined on a partition 
of Cl. Let Vh denote a partitioning of Cl into Nk = Nx(Vh) subdomains I\ with boundaries 
dK such that 

(i) N K (V h ) < oo 

(ii) H = U{A' : K € Vh] 

(iii) For any pair of elements K, L € Vh such that K ^ A, K fl L = 0 

(iv) K are Lipschitzian domains with piecewise smooth boundaries 

(v) dK- = {x € dK\ f3 ■ n/c < 0} and dK+ = d K \ dK- 

(vi) T h _ — U^:=i dK fl r_ coincides with T_ for every h > 0 

(vii) Tkl = dK fl dL is an entire edge of both K and L 

Let V(K) = {iME L 2 (K)\ v 0 e L 2 {I<)}, then V(V h ) = n^i V(I<). Note that a function 
v € V{Vh) need not be continuous across element interfaces. We use the following notations 
concerning functions v,w (E V{Vh) '• 

v mt h = u| A .(x), x e dK 
v exlK _ v \ l ( x ), xe dKf]dL 
v ± = limu(x±e/3) 
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(v,w)-y = / vw\j3 ■ n e | ds 7 C 9/\ 


((»>)i = (»>»)- 


Ml 


Nk Nk , 

= E IMI 2 a-= £ L v2dx 

K=1 K = 1 


The problem corresponding to (l)-(2) defined on any partition Vh is as follows: 


Find u(x) £ V'X'P/,) such that for every K € Vh 


up -{- au = f in K 

u ml h /3 ■ uk = u extft ^-nK Vx € d/F \ dfi > PI 
ii int A /3 • n^- = g j3 • nA Vx (E 5 K fl T_ 

Let P2 denote the following variational boundary value problem for any partition Vh'- 


Find u(x) £ V(Vh) such that 
B(u,v) = C(v) V v £ V(Vh) 4 


P2 


where 

N k 

B(u,v) = ^2 {(up + au,v) K + (u + - u-,v+) dK _\r_ + (u,v) dK - nr_} (4) 

K=l 

n k 

£{ v ) = £{(./»* + (^v)aiifjnr-} ( 5 ) 

A'=l 

Note that with the definition of V('Pfc) we can apply Green’s formula to (4) to obtain 

Nk 

B(u,v) = J2{(u,-vp + av) K + {u ,v - v + ) dIi _\r_ + {u, v) dK+ nr+ } ( 6 ) 

K=\ 
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Lemma 1 Let the bilinear form B( •, •) be defined by (f). Then there exist positive constants 
a and M such that 


B(v,v) > allHHl VvEV(Vh) 

B(v,w) < M|||u|||/j IIIHII/J Vv,w6V(V h ) 


where 


Ill-life = E{IMfe + ((- + — ■ + <W)^n»} 

K = 1 

Nk 

Ill-life = E{ll-olli + ll-lfe + «- + )&_\r_ 

A'=l 

+ (( v ))dh'-\ r_ + {( V ))9/Cnafi} 


Proof: (i) From (4) we have 

B{v,v ) = i( v 0 + qv i v )k + ( v+ ~ v~,v + ) 9K _\ r _ + (v,v) dK _ nr _} 

Kev h 

Applying Green’s formula to the first term: 

{vp, v) K = l ~ J K (v 2 )pd x =\j dK v 2 P ■ n A ds 

= 2^ V+ ^ dK ~ ^ ° ^ dK + ^ 

and noting that 

N k N k 

E «-‘»k = E {«--»J K _\r_ + ((-■))«* nr t ) 

A'=l A'=l 

yields 

Nk i i 

B(v , u) = Y, {( GU ’ v )k + 9 (( v+ - v ))Ik.\ r_ + 9 « u ))L + nr + } 
K=l Z Z 


(7) 

( 8 ) 


(9) 


( 10 ) 
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from which the first inequality follows. 

(ii) Adding the definitions of £?(•,•) in (4) and (6) yields 

N h - 

2B(v,w) = {( v 0’ w )k ~ ( v , w p)k + 2(a«,u?) K 

K=l 

+ (v + - v~ ,w + ) dK _\r_ + (v~,w~ - w + ) dK _\ r_ + (u,w)aA'nan} 


Applying the triangle inequality to the jump terms and the Holder’s inequality to the inte- 
grals and resulting sums yields 


, rr K 

B(v,w) < imax(l,2||a|U n ){S[IMIt+2|Hlt + ((« + ))L_\r- 

K=l 
2 


+ 3 (( v })gA'_\r_ + (( u ))aifnan ]} 2 


x { £ lIMli + 2|Hli + 3 (( w+ ))L_\ r. 

K = 1 

+ (( w ))aA"_\r_ + (( w ))9A'nan]} 2 


from which the second inequality follows. 

Remark: We note that it is sufficient to take the constants in the bounds of the Lemma to 
be the numbers 


a = minfmina(x), -) 
v xgn v ; 2 

(11) 

3 

M = -max(l, 21101100,0) 

(12) 


I 


4 Discontinuous Galerkin Approximation 

Approximate solutions to P2 are sought in a finite dimensional subspace of V(Vh) which we 
denote by V p (Vh) and define precisely below. The discontinuous Galerkin approximation is 


7 



obtained by replacing u,v <E V('P/ l ) by 6 V p (Vh) as follows: 


Find u p h 6 V p (Vh ) such that 

5K,0 = £«) (13) 

where B(u p ,v p ) is given in (4) and C(v p ) is given in (5). 

4.1 The Finite Dimensional Space V p (Vh ) 

Let the elements A' € be quadrilateral elements which are affine maps of a master element 
K = [—1,1] x [—1,1], i.e., I\ = Fk(I\) as illustrated in Fig. 1. Let h k = diam(A'), S be 
a sphere contained in K, and pa = sup{diam(S’)}. For the analysis we assume that Vh 
belongs to a family T of quasiuniform refinements, that is for every K (E Vh, there exist 
positive constants a and t independent of h — max^ep h hh such that 

< t and < a (14) 

hh' PK 

The finite dimensional space V v {Vh) C V(Vh) is defined as follows 

V p (Vh) ={ve L\ SI) : v\ K 0 F k = v K € Q Ph (K)} 

where Q Pk (K) is the space of tensor products of complete polynomials of degree < p K 
defined on the master element K. We use the notation that v K £ Q Vh (K) to imply that 
v K € Q Pk (K). The basis for Q Ph (K) is formed by tensor products of one-dimensional 
Legendre polynomials. Note that in general the elements of V p (Vh) are discontinuous across 
element interfaces and that the degree of the polynomial approximation may vary from one 
element to the next. 

In proving a priori error estimates for solutions of (13), we will need the following basic 
approximation properties of functions belonging to V p (Vh)- Since functions in V p (Vh) are 
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discontinuous at element interfaces, we are primarily concerned with polynomial approxima- 
tions on a single element and its boundaries. 

Lemma 2 Let I\ € Ft? be an affine map of a master element K = [— 1 , 1 ] x [— 1 , 1 ], that 
is K = Fk(K). Let 7 denote any edge of dK which is an affine map of a master edge 
7 = [—1,1], Let wk be a polynomial of degree Pk defined on the master element. Let 
wk = wr 0 Fk denote the image of wk under the transformation Fk- Then (3 ■ V% 
satisfies the following: 

II/J-Vu.kIIk < C^IKIIk ( 15 ) 

UK 


J (l 3 • Vw K ) 3 

1/3- 

n-y| ds < 

CjF. « w K »* 

( 16 ) 

where the constants C are independent 

of h K ,p K , 

and wk- 


Proof: For polynomials of degree 

PK 

on the master element we have that (see 

Dorr [ 7 ]) 


< 

IK-IU 

< Cpk\\w k \\k 

( 17 ) 

| WR |s,7 

< 


5 : ^Pk\\wk\U 

( 18 ) 

where the constants C > 0 depends 

on 

s , but not 

on p K or WK' 



For affine mappings Fk, a standard scaling argument (see Ciarlet [8]) yields that for s > 0 
an integer, there exist constants C > 0 such that 



< 

Ch l K 3 | l f’A'ls,A' 

( 19 ) 

|wjf|.,7 

< 

Ch^ K |wa'|*,7 

(20) 

I^A'I s,k 

< 

Ch g K l \w K \ s , K 

(21) 


< 

ChK>\w K \,„ 

(22) 


where C depend on s, a, and r (see ( 14 )), but not on h K , p K , or wk . 
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The first estimate (15) follows by combining (9), (7), and (21). The second estimate (6) 
follows from (20), (18), and (22). B 

We also have the following result from Babuska and Suri [9] concerning the polynomial 
approximation of functions on a. single element. 


Lemma 3 (Babuska and Suri [9]) Let K £ Vh, 7 denote any edge of dK , and u £ 
H S (K). Then there exists a constant C = C(s,r, a) (see (14)) independent of u,p k, and 
h K , and a sequence z v h £ Q Pk ( I<),p K = 1,2,... such that for every 0 < r < p K 


u 


u 


Pk 

s > 0 

(23) 


Pk 2 

1 

^2 

(24) 


where v = min(pA' + 1 , 3 ). 

4.2 A Priori Error Estimate 


The discontinuous Galerkin method (13) was first analyzed by Lesaint and Raviart [1] for a 
given fixed value of p K , i.e. for the case in which p K = p for every element K £ Vh- The 
error in a solution Uh to (13) approximating an exact solution u £ H 3 (Cl) to P2 was shown 
to be 

Il u ~ < C^ -1 ||«lkn 

This estimate is not optimal in the sense of interpolation error estimates and was improved 
by Johnson and Pitkaranta [2]. Using a mesh-dependent norm they showed that 

|||tl-tlfc||| M <C&-2||«|U 


where 


\h,p 


ll« - w ^lll 


2 

h,P 
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- X ^/dl e /?lli 2 (K) + I l e l li 2 (fi) + X (( e+ ~ e ))h<-\ r_ + (( e ))an 

A'eP h Ker h 

While this estimate is not optimal in the sense of interpolation error estimates for j|e||^ 2 (Q), 
it is optima] with respect to \JFk\\z()\\l' 2 (K) and {( e+ — e- ))a/c_\r_- 

We shall derive estimates similar to Johnson and Pitkaranta [2]. Taking into account 
that p K is not constant, we shall use the following mesh-dependent norm 

IIMIIa p /= X {^rlMI K + I Ml* + «^ + -v~))l K _ v _ + (25) 

Kev h Pk 

The presence of the polynomial degree p K in this norm deserves comment. For each element 
I\ E Vh we assign to K a positive integer X K which serves as a weighting factor tuned 
to allow optimality in some sense of the error estimate in a mesh-dependent norm. Later, 
the choice A /c = p K will prove to be the proper one since it allows for the development of 
quasi-optimal hp-estimates. 

We first prove that the bilinear form in (13) satisfies an inf-sup condition with respect to 
this norm on the space V p (Vh ) and therefore that the solution depends continuously on the 
data. 

Lemma 4 For every v p h E V p (Vh), there exists a w p E V p (Vh) such that 

B( v hi w D ^ C B \M\\\l P ,p (26) 

IIKIII/ip.0 < £||I u EIJUp>0 (27) 

where the positive constants C B and C are independent of h K , p K , and v%. Moreover, the 
solution u p h to (13) satisfies the following global stability estimate: 

IIKIIUp,/? < C'dl/UnT << g >>r_) 
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Proof: Define the restriction of w p h 6 V p {Vh) to an element K € Vh as 


P K 


( 28 ) 


where 8 £ (0, 1] is defined later in the proof. Dropping the h,p, and K scripts for ease in 
notation, we have 


Bk(v,w) = j {vp av){v + 8 -j-vp) dx 
Jk Pfc 

+ / (v + - v~)(v + + 8 ^-vp + )\(3 - n/,1 ds 

JdK-\ r_ pi- 


+ f v + (v + + 8 ^-vp + )\f 3 ■ tik\ ds 

JdK-nr- p K 


> a 0 \\v\\ 2 K + 6 ^-\\vp\\ 2 K + I vvp dx 
PK Jk 


2 

y dK-\ r_ 


+ 6 -- f vvp dx+ « v + >>g R -_nr_ + << v+ >> \ 

Pk Jk 

~ L.\ r_ V * V ~ V ' n,fl ds + s! jstL-V- {v * ~ ' nA '' * 

+ sh £.L-n 


Noting that 


/ m, dx=lf (v-y\D • n K | is - i / (<> + ) 2 l0 • n K | <fe 

«/ A ^ */ ^ JdK— 


(29) 


and that from Lemma 2 


I j K wp dx\ < c^WvWl 

| f v + vp + \(3 ■ n/v’l ds\ < ci << u + >> 
A "‘K 
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we have 


Bk(v,w) > (a 0 ~ q£)IMIa’ + ^^IMIa 

Pk 


+ 2 <<: v+ >> dK-\ r_ +(^ — C2 ^) << v >>2 di<-nr- 

- [ v + v~\/3 • njf| ds + 6^y [ {v + - v~)v 0 + \fi ■ n K \ 

JdK-\ r_ pi- JdK-\r- 


ds 


Using the Schwarz inequality and the previous inequalities, one can show that 

\^~r I ( v+ “ v ~) v i3 + \P ' n A'l ds \ ^ ~^(« u + »k_\r_ + << v ~ »k-_\ r_) 

P K JdK-\r _ 1 

Now summing over all the elements K (E Vh and realizing that 
{ 2 « v ~ »Ik-\ r_ — « v >>2 dK-\ r_ 

Kev h 

+ (g 2”^) << v+ >>2 dK-\ r_ ~ J gK u+t; ‘ nA 'l ds } 

> l « v »l + + min(l, i - ^y<5) « n + - »k_ 


\r- 


results in 






«KX) > («o-c«)IMlS + < £ -flMI* 

A'G'Ph ™ 


+ (5 “<**)«»»£_ +\ «n»r + 

+ min(l ^ ~ -7T 6 ) H « n + - v~ »k'_\r_ 
z z A'en 

Choosing (5 = min(|, g^) yields the first inequality. 

The second inequality easily follows from the definition of w v h and Lemma 2. The third 
inequality follows by combining (26) with (13) and (5): 

C B \ \\<\\\l p , 0 = B(ulvl) = C(vl) 
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where v£| K = u p h \ K + <5/3 • Vu£| fv . and applying Holder’s inequality to C(v p h ) defined in (5). I 
We now have all the preliminary results needed to prove an a priori error estimate for 
an hp - version of the discontinuous Galerkin method. 


Theorem 1 Let u € H T (Ll) be a solution to P2 and let u v h be a solution to (13). Then there 
exists a positive constant C independent of h K , p K , and u such that the error, e — u — u p h , 
satisfies the following estimate 


\e\\h P ,0 < C \ 5Z 

,K€V h 


hp '- 1 

n2r — 2 


max 



(30) 


where v K — min(p K + l,r). 


Proof: Let n v h u € V p (Vk) be an /ip-approximation of u that satisfies the estimates in Lemma 
3 and write 

e = u — u p = u — n£u + — u v h 

which implies that 


WHU < 2(|||u - njuiHijy, + IIK - njuiii^) 


(31) 


Subtracting (13) from P2 yields the orthogonality condition that 


= B(u - nju, »;) - B(ui - nj«, t>Z) = o 


v«2 € V r (V>) 


Combining this with Lemma 4 yields that 


cjik - n;u|||J^ < bk - nj»,«) = b(« - ns»,«) 


( 32 ) 


for a particular choice of v € V v (fPh )• To simplify the notation, let 


rj = u — U p h u and w = u p h — II p h u 
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and recall from Lemma 4 that the particular choice of v € V p (Vh ) for which (32) holds also 
satisfies the estimate 


IIHIk, < Clllu.HI*,,, ( 33 ) 

Next we seek to bound from above B(rj, v) on the right hand side of (32). Using the definition 
in (6), we have 

B{v,v) = ^2 {(r],-V 0 + av) K + - v + ) aK _\ r _ + {r),v) dK+nr+ } 

KzP h 


^ v~* Pk || || y/fiK 

< L -Tr^lmh' ■ ~zr~ 

KeP h V h ' K 


PK 


Mk- + 


I oo , I\ 


K*A- 


+ {(*? ))aA_\r-((^ + “ ” ))aA'_\r_ + (( T l))dK + nr + ((^))dK + nr + } 


Using Holder’s inequality for sums results in 
B(t],v) < max(l, ||a||oo 4 j) 


N 


£ p hw\l 


KeP h ' 


E^INIi + ./Elhlll./ElWR 


"ill** Z-f „2 

K \i<ev h Pk 


Kev h 


KePh 


+ ,/ E ((*? ))dK-\r-J ((^ + ^ )) 9 A-_\r. 


KePh 


KePh 


+ ./ (( 7 /)) 9 A' + nr + . I ((^))aA' + nr + 

V KePh V KeVh 

Applying Holder’s inequality for sums again yields 


B(v,v)<C\J 2 


Pk 

K L h K 


Combining (35), (33), and (32) results in 


|<-HHIU^<c £ 


Pali 

&A' 


^IIa + Ml + i(v ))aA'+\r + + (( P))dK+nr+ 


Using the estimates (23)-(24) in Lemma 3, we have 


h v K 

K < C^-WuWrJC 

Pk 


(34) 


(II* + HA + <fo - »S*r+\r* + <W>Jjf t nr t ]} ! lll*lll»w» (35) 


( 36 ) 
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"A- 1 
A 


h,7 


\M, <C((3- n K )M\ hK < C^-IMU- 

P K 


«r» 8Kl \ r, < 

Pa- 2 

where v K = min(p K -f l,r). 

Substituting the above estimates involving r) into (36) yields 


IK-iMM < c{ E 

,KZV h 


h)y~ x h 2 ^ h]- K ~ l 

Pk ~ 2 + Pa + Pa " 1 ' " U| 


r,K 


(37) 


Recalling (31), we see that all that remains is to bound |||u — Il£u| \\h P ,p which follows from 
Lemma 3 and (25): 


lll«-niu||| w <c{ E 

[Kev h 


,2v k ,2u k -1 ,2 v k 1 
n K , 11 K , 11 K I II i|2 


Combining (31), (37), and (38) completes the proof. 

Remarks: 


(38) 

■ 


2* k -1 

(i) For < 1, the estimate becomes |||e|||* Pt p < C{Y,Kev h ll u Hr, k ) 5 

(ii) For p K =constant, the a priori error estimate reduces to the one derived by 
Johnson and Pitkaranta [2]. 

(iii) The error estimate reveals that the discontinuous Galerkin method provides some 
natural control in the /3-derivatives of the approximate solution. The factor 
means, however, that this control decreases as ^ — ♦ 0. 

5 Numerical Examples 

We verify the estimate in (30) with two examples where j 3 = (0.8,0.6) T , a(x) = 1.0, and 
ft = (-1,1) x (-1,1). 
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Figure 2: Quadrilateral element mesh used for quasiuniform refinements in Example 1 

5.1 Example 1 

In the first example, the source term /(x) was chosen so that the exact solution to (1) is the 
C°°(Q) function 

u(x) = 1 + sin ^(1 + x)(l + y) 2 ^j 

with an inflow boundary condition of g = 1. The error in the solution obtained with 
varying h and p is listed in Table 1 for uniform refinements of a mesh consisting of square 
elements. The error in the solution obtained for quasiuniform refinement of a mesh consisting 
of quadrilateral elements (see Fig. 2) is listed in Table 2. 

To verify the estimate (30), we first consider the case when p K is fixed and h K is varied. 
According to (30), we should get |||e|||A P ,/j < Ch^ c+a ||u|| r> n. This is verified in Fig. 3 
where |||e|||^ is shown as a function of h K . On the log-log scale, the slope of the lines 
corresponding to a fixed value of p K is indeed p K + | for both the uniform and quasiuniform 
meshes. Next we consider the case when h K is fixed and p K . is varied. In this case, the estimate 
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Mesh 

- log h 

-log || u - uaIIUp,/? 

p = 1 



I 

2x2 

0.000 


[MU 



4x4 

0.301 

0.5552 

[ 



8x8 

0.602 

0.9692 




16 x 16 

0.903 

1.4003 

3.1163 

4.5834 

6.3047 

32 x 32 

1.204 

1.8412 

3.8574 




Table 1: Example 1 _ Error using uniform hp meshes 


Mesh 

log h 

-log |||u - Uh 1 1 |hpj3 

P = 1 

p = 2 

wmmm 

ummsm 

p = 4 

2x2 

-0.2116 


0.8586 

1.7402 

2.2831 

4x4 

0.0689 

0.5153 

1.5930 

2.5395 

3.4998 

8x8 

0.347 



3.5814 

4.9723 

16 x 16 

0.641 

1.3913 

3.0955 

4.6208 

6.3196 

32 x 32 

0.938 

1.8129 

3.7870 




Table 2: Example 1 _ Error using nonuniform h and uniform p meshes 
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Figure 3: Example 1- Rate of convergence of error for fixed p 


(30) reduces to |||e|||/i P>j g < Cp K ~ r + l ||u[| ri n- Since u G C 00 ^), we should expect exponential 
rates of convergence. This is confirmed in Fig. 4 where the curves corresponding to |||c|||^ Pt /j 
as a function of p K have a slope on the log- log scale which increases as p K increases. These 
results are combined in Fig. 5 where 1 1 |e| | is shown as a function of the total number 
of unknowns in the solution. The solid lines represent /i-refinements for a fixed p and the 
dashed lines represent p-enrichment for a fixed h . Clearly for smooth solutions, higher-order 
accuracy is achieved for the same number of unknowns by using higher-order elements, 

5,2 Example 2 

In this example the source term /(x) was chosen so that g = 1 and the exact solution to (1) 
is the C°°(Ct) function 

U(x) = 1 + ^(1 +*)(! - *)(1 +y)(l - J/)tan -1 (o;(£ — £o)) 
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p 


Figure 4: Example 1- Rate of convergence of error for fixed h 



Figure 5: Example 1- Rate of convergence of error with number of unknowns 
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Mesh 

-log h 

-log || u - u h \\\ hPi! 3 

p = 1 

P = 2 

II 

CO 

p — 4 

2x2 

0.000 



1.7957 

2.1453 

4x4 

0.301 


1.9235 

2.5461 

3.0656 

8x8 

0.602 

1.5571 

2.6252 

3.4634 

4.2565 

16 x 16 

0.903 

2.0017 

3.3471 

4.4822 

5.5530 

32 x 32 

1.204 

2.4512 

4.0925 




Table 3: Example 2 _ Error using uniform hp meshes 


where 

£ = 2 + Z j~ y , o = 10, and £ 0 = 0.6 

2y2 

The error in the discontinuous Galerkin solution obtained on uniform meshes with various 
values of h and p are listed in Table 3. The error in the solution for uniform ^-refinements 
with fixed p is shown in Fig. 6. The error for uniform p-enrichments with fixed h is shown 
in Fig. 7. The error as a function of the total number of unknowns is shown in Fig. 8. The 
rate of convergence of the error given in Theorem 1 is verified for this example. 

6 Concluding Remarks 

The discontinuous Galerkin method can be viewed as an elementwise application of the 
standard Galerkin method in which jumps on element boundaries are admitted naturally in 
the formulation. The method can also be viewed as a natural higher-order extension of finite 
volume methods, except that the coefficients in the higher-order polynomial representation 
of the solution are obtained by solving the conservation law and not by some post-processing 
of solution mean values. 

In this paper, we have derived an a priori error estimate for an hp - version of the discon- 
tinuous Galerkin method. These estimates are derived in a mesh dependent norm in which 
the coefficients depend upon the local mesh size h K and a number p K which can be identified 
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Figure 6: Example 2- Rate of convergence of error for fixed p 






Figure 8: Example 2- Rate of convergence of error with number of unknowns 

with the spectral order of the local approximation over each element. The estimate is general 
in the sense that it is valid for general quadrilateral elements of any degree p K . Moreover, 
the estimate is valid for meshes in which the size and the spectral order of the elements vary 
throughout the mesh. 
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